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Abstract: Optical whispering gallery mode (WGM) resonators have been 
very attracting platforms for versatile Kerr frequency comb generations. 
We report a systematic study on the material dispersion of various optical 
materials that are capable of supporting quality factors above 10 9 . Using an 
analytical approximation of WGM resonant frequencies in disk resonators, 
we investigate the effect of the geometry and transverse mode order on 
the total group-velocity dispersion ( GVD ). We demonstrate that the major 
radii and the radial mode indices play an important role in tailoring the 
GVD of WGM resonators. In particular, our study shows that in WGM 
disk-resonators, the polar families of modes have very similar GVD , while 
the radial families of modes feature dispersion values that can differ by up 
to several orders of magnitude. The effect of these giant dispersion shifts 
are experimentally evidenced in Kerr comb generation with magnesium 
fluoride. From a more general perspective, this critical feature enables 
to push the zero-dispersion wavelength of fluorite crystals towards the 
mid-infrared (mid-IR) range, thereby allowing for efficient Kerr comb 
generation in that spectral range. We show that barium fluoride is the most 
interesting crystal in this regard, due to its zero dispersion wavelength 
( ZDW ) at 1.93 /im and an optimal dispersion profile in the mid-IR regime. 
We expect our results to facilitate the design of different platforms for Kerr 
frequency comb generations in both telecommunication and mid-IR spectral 
ranges. 
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1. Introduction 

Optical frequency combs have been very useful scientific and technological tools in a wide 
range of applications such as optical clocks, precision spectroscopy, low phase noise microwave 
generation and optical communications. The general approach for producing frequency combs 
is based on mode-locking in lasers. An alternative approach has been recently demonstrated in a 
silica monolithic toroidal microresonator (IJ, where cavity enhanced four-wave mixing (FWM) 
occurs. As a result, this type of comb using the third-order susceptibility in resonators is usually 
referred to as Kerr optical frequency comb. The emergence of such combs provides a compact 
and energy efficient solution. Therefore, extensive amount of works in both theories and exper¬ 
iments have been devoted to this field. For example, the understanding of the comb formation is 
now very elaborate, and involves both coupled-mode and spatiotemporal models (2}j7]]; various 
applications have also been experimentally explored ranging from low phase noise microwave 
oscillations & ultra-short pulse generation (10|[TT), and optical communications CD- 

Compared with Kerr comb generation platforms based on chip-scale waveguide racetrack 
resonators CD , the mechanically polished disk resonators have a much wider range of material 
choices and can feature comb generation with very low pump power in the sub-mW level 0. 
which results from the ultra-high quality ( Q ) factors and small mode volumes in these WGM 
resonators. Although Kerr frequency combs in fluoride WGM resonators have been demon¬ 
strated in the telecommunication window and recently at 2.5 /im in the mid-IR regime CD’ 
there is still not a systematic study on the optimal choice of materials and the sizes of res¬ 
onators. Moreover, anomalous GVD is generally desirable for comb generation and the span¬ 
ning of Kerr combs is also strongly dependent on the dispersion profile. Recently, an analytical 
approximation for the dispersion of WGMs in the disk geometry has been derived by two dif¬ 
ferent groups, showing a good agreement with finite-element simulations and high-precision 
free spectral range (FSR) measurements Q5j[l6j. It enables an accurate study of the impact of 
geometry on the GVD. On the other hand, barium fluoride (BaF 2 ) with a good transmission 
range covering the whole mid-IR window has been recently shown to be capable of supporting 
Q factors above 10 9 CD , and this result broadens the choice of crystalline candidates for Kerr 
comb generation. 

In this article, we report the study of material dispersion for different IR and mid-IR optical 
materials using Sellmeier equations. We point out that fluoride crystalline materials generally 
have significantly flatter material dispersion profile than other materials. A systematic and the¬ 
oretical investigation of the impact of resonator geometry and WGM indices on the GVD is 
carried out, showing that the major radii and the radial WGM indices are key parameters to 
consider for GVD tailoring in WGM resonators. The impact of geometry on GVD is also ex¬ 
perimentally evidenced in an experiment showing how primary Kerr combs with significantly 
different multiplicity can be generated with the same MgF 2 resonator, but pumped in WGMs 
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Fig. 1. (a) Geometry of the disk resonator, (b) Schematic of Kerr frequency comb genera¬ 
tions in a fiber taper coupled resonator setup. WGMR: WGM resonator, vq is the primary 
parametric gain peak spacing from the pump frequency vq. 


with different radial orders. We also investigate theoretically MgF 2 and BaF 2 platforms for 
potential broadband comb generations in the telecom and mid-IR regimes. 


2. In search of the optimal materials for Kerr comb generation 

As far as Kerr frequency comb generation is concerned, the primary step has been the search 
of the optimal materials. Although the on-chip ring resonator as comb generation platforms 
is much more compact, mechanically polished high Q WGM crystalline resonators are still of 
great interest, due to the flexible choice of materials and ultra-high Q factors leading to power 
efficient operation. 

In these axi-symmetric disk resonators, the rim part can be treated as a toroidal boundary 
with the main radius R and the curvature radius r, as illustrated in Fig.[l](a). Figure |T](b) gives 
a typical schematic drawing of Kerr comb generation. A single frequency continuous-wave 
(cw) pump laser can be coupled into the WGMR through a fiber taper or a prism (T8). Under 
certain conditions that have been analyzed in refs. (2j[3]|7), the built-up intracavity light leads to 
the modulation instability (MI) gain with offset vg from the pump frequency vo. The parametric 
FWM then creates the first two sidebands where WGMs with single or multiple free spectral 
range ( FSR ) relative to the pump mode meet the MI gain. Subsequent parametric process cre¬ 
ates a comb with multiple-FSR spacing, generally referred to as the primary comb. Depending 
on the dispersion and pump condition, broadband Kerr combs with single-FSR spacing and 
thousands of frequency lines can be further produced. As Kerr comb arises from MI gain, it 
thus prefers an anomalous dispersion profile similar to MI gain in optical fibers. In addition, 
the MI theory in Kerr comb generation predicts that the offset vq of the parametric gain peak 
frequency scales with the dispersion following vg 1 /VGVD [2,3,7 ]. Hence, weak dispersion 
is generally favorable for broadband Kerr comb generation. It should be noted that the primary 
Kerr comb shown in Fig.|T](b) also yields rolls (or Turing patterns) in the spatial domain f7p9| . 

At this date, Kerr optical frequency combs have been reported on WGM resonators made of 
fused silica, magnesium fluoride, calcium fluoride (CaF 2 ) and sapphire |lj] 14 20-24), which 
are ah capable of featuring Q factors above 10 9 . Recently, such high Q factors have also been 
demonstrated on barium fluoride CD- We thereby take into account this material with others 
for the investigation of the material dispersion ( GVDm ) using their Sellmeier equations (25) . 
The GVDm is usually characterized by group delay dispersion (in ps/km/nm): 
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Fig. 2. Comparison of material dispersion for different optical materials: fused silica 
(Si02), sapphire, MgF 2 , CaF 2 and BaF 2 . Dashed lines designate the corresponding ZDW s : 
1.27 gm, 1.31 gm 1.34 jUm, 1.55 gm and 1.93 gm, respectively. The vertical blue and 
brown zones highlight the telecommunication C-band and mid-IR windows. Note that the 
GVDm profile of MgF 2 investigated here is derived for the ordinary ray. The extraordinary 
one has very similar profile with a slightly different ZDW at 1.57 gm. 


where A is the wavelength, c is the speed of light in vacuum, n(?i) is the refractive index. 

Figure [2] displays the corresponding material GVDm as a function of wavelength from 0.8 to 
6 gm. It should be noted that the dispersion is anomalous when GVD is greater than zero. A 
broadband comb generation will require a small GVD value in this regime. On the other hand, a 
GVD that is too large will also make Kerr comb generation difficult, as the corresponding para¬ 
metric gain could be too close to the pump frequency and do not meet any WGMs (the phase 
matching condition will therefore never be fulfilled). One can also see that fluoride materials 
have flatter dispersion profiles than the others. The corresponding ZDW of these materials are 
listed in Table [T] Also shown are some optical characteristics of these materials (26). 

One can note that fused silica, quartz and MgF 2 have anomalous GVDm around the tele¬ 
com wavelength, while all the fluoride crystals are transparent well into the mid-IR regime and 
also feature anomalous GVDm . Among them, BaF 2 is characterized by a flat high transmis- 


Table 1. Characteristics of bulk fused silica, Sapphire, MgF 2 , CaF 2 , BaF 2 at T ^ 300 K for 
Kerr comb generation a 


Name 

Transmission 

(Mm) 

Ao 

(Aim) 

Refractive 
index at Ao 

(dL/L)/dT 

(ur 6 /K) 

dn/dT (10~ 6 /K) 

Fused silica 

0.2-4.0 

1.27 

1.4473 

0.55 

9.6 

Sapphire 

0.19-5.2 

1.313 

1.7502 

6.65(a)/7.15(c) 

13.6(o)/14.7(e) 

MgF 2 

0.13-7.7 

1.34 

1.3717 

9.4(a)/13.6(c) 

0.9(o)/0.3(e) 

CaF 2 

0.12-10 

1.55 

1.4260 

18.9 

-11.5 

BaF2 

0.14-13 

1.93 

1.4648 

18.4 

-16.2 


a ( dL/L)/dT along the a and c axes are denoted, o: ordinary ray; e: extroordinary ray. 
dn/dT for fluoride crystals and sapphire are given at 1. \5\im and 0.589 gm, respectively. 
dn/dT for fused silica at 1.13 jlm is taken from manufacturing specifications. 























sion covering the whole mid-IR regime (2.5 — 10 /im). Moreover, the reported 2 nm surface 
roughness in BaF 2 Q7), also confirms its potential for reaching ultra-high Q factors that are 
only limited by the material absorption across the whole mid-IR regime. Also shown in the 
table are their thermal expansion coefficients ( dL/L)/dT and thermo-optic coefficients dn/dT. 
These parameters are particularly important for WGM applications. The self-thermal locking 
method can be used for WGM pumping for microlasers (27) and combs Q- However, it is 
only applicable for the material that possesses the same signs of these two parameters, such as 
MgF 2 , fused silica and Sapphire. On the other hand, the others could benefit from self-thermal 
compensation like athermal glasses. Hence, they are less sensitive to the environment tempera¬ 
ture fluctuation, due to the fact that the frequency variation dv of a WGM is determined by 
(dv/v)/dT « (dn/n)/dT — ( dL/L)/dT [28]. Considering another useful locking technique, 
the self-injection locking method using Rayleigh scattering induced feedback is applicable in 
both materials (9j[29|. 

3. Geometry dispersion effects of WGMs in disk resonators 

It is known that the geometry dispersion can have a strong impact in the total dispersion, such 
as that of the silicon channel waveguide (30) . Similarly, it also plays an important role in res¬ 
onator structures. As light in WGM resonators is confined under the equatorial boundary of the 
disk, the size and the curvature of the rim are expected to affect the total dispersion value. A 
recent study also shows that the dispersion of a spherical resonator can be tailored by a coating 
layer 0 Here, we focus on monolithic resonators without any coating to preserve their high 
Q factors that are only limited by the absorption of the host material. 

A simple analytical approximation of the eigenfrequencies of WGMs derived by two 
groups (T5j[l6} will be used in this work. The WGM resonant wavelength in a toroid can be 
expressed as follows: 
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where n(X) is determined by the Sellmeier equation of the material, m, q,p are the mode num¬ 
bers in the azimuthal, radial and polar directions, a q is the negative qth root of the Airy function. 
The index p is similar to that of WGMs in spheres and is related to two momentum numbers / 
and m by p = / — m, indicating the p + 1 antinodes of the field distribution in the polar direction 
which can be identified through excitation mapping (32). Other methods from spectroscopy 
and emission patterns for identifying WGMs have also been demonstrated [ 33p4| . In a typical 
mm-size disk resonator, m usually ranges from 10 3 to 10 5 for optical resonances. 

It should be noted that Kerr combs have been generated using two different sets of WGMs. 
The first one is the family of WGMs with spacing of one FSR, vfsr = v m — v m _i and the other 
one is a group of polar WGMs with v p — v p -\ spacing |21 ], where v = c/X is the resonant fre¬ 
quency. Here, we focus on the Kerr combs generated using the longitudinal FSR, and which can 
feature hundreds to thousands of comb lines, interesting for many spectroscopy applications. 
We therefore analyze the total dispersion of these longitudinal FSR spaced WGMs. The disper¬ 
sion of WGMs underlying Kerr comb generation is characterized by the variation in the FSR, 
Avfsr = (v m +1 — v m ) — (v m — v m —i ). This value can be derived by numerically solving Eq. [2] 
with the underlying Sellmeier equations. To convert Avfsr into the widely used dispersion pa¬ 
rameter GVD in units of ps/nm/km, we applied the following approximated expression (35): 
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Fig. 3. (a) Calculated geometry dispersion influence on the total GVD for a MgF 2 disk with 
different R and r = 0.1 mm. Dashed line: material GVDm curve, (b) Calculated zero total 
dispersion wavelengths Ao for R from 0.5 mm to 5 mm. Inset: distribution of the electric 
field amplitude for a fundamental WGMs with R = 1 mm and r = 0.1 mm. Note: Ao is the 
resonant wavelength that has its total dispersion closest to 0. 


Here, we will study the impact of geometry dispersion by assuming WGM resonators made 
of MgF 2 , which has a small GVD over a large range when compared with fused silica and 
quartz. We first consider the case of the fundamental WGM family with q = 1 and p = 0 and 
discuss the size impact of R and r. As one can see in Eq. ([2]), the second term on the right side is 
much larger than the third term due to the large value of m in mm-size resonators. We thus study 
only the impact of the main radius R on the total dispersion, because the effect of r is similar to 
the effect of the polar number p that will be introduced later. Figure [3] (a) shows the calculated 
total GVD for different radii R values from 0.5 mm to 5 mm around the zero dispersion regime. 
The dashed line presents the material GVDm . 

One can note that smaller R shifts the zero GVD further into the longer wavelengths, which 
results from a normal geometry GVD adding on to the material GVDm . The total GVD profile 
compared to the material one is eventually pulled down. Similar effects have also been reported 
[ 23,36 ]. The wavelengths (Ao) of the WGMs closest to zero GVD as a function of the radii R is 
shown in Fig.[3](b). The inset provides a calculated electric field distribution of a fundamental 
WGM resonating at Ao in a MgF 2 disk with R = 1.0 mm and r = 0.1 mm using an analytical 
expression given in G3 From this figure, one can conclude that small anomalous GVD can 
be engineered by carefully choosing the major radius of the resonator. This is a key step for 
designing a Kerr frequency comb setup. For example, the GVD for the fundamental modes in a 
R = 0.5 mm MgF 2 resonator at 1.55 /im is reduced to 0.72 ps/km/nm, which is about one tenth 
of the original material GVD (7.0 ps/km/nm). Careful engineering of this parameter, the GVD 
reduction can be even larger. 

The ability of the GVD tailoring in a WGM resonator after fabrication is also critical. 
Thereby, we have also investigated the effect of different transverse mode orders on the to¬ 
tal GVD. Figure]?] shows the corresponding calculation for a MgF 2 resonator with R = 1.0 mm 
and r = 0.1 mm. The black dashed lines marks the material dispersion itself. The upper images 
are the corresponding calculated field distribution patterns. The impact of WGMs with different 
q indices on the total GVD is presented in Fig.[?](a). On the other hand, the mode index p plays 
a less interesting role in tailoring the total GVD as shown in Fig. [4] (b). The inset is a zoom in 
of the GVD profile around the ZDW. The red dashed line highlights the telecom wavelength of 
1.55 /im. The fundamental radial mode with q= 1 possesses a GVD of 3.0 ps/km/nm, which is 
much larger than the previous value for the fundamental mode in a resonator with R = 0.5 mm. 
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Fig. 4. Calculated total GVD for different transverse WGMs in a MgF 2 disk with R = 1 mm 
and r = 0.1 mm. Top: Electric field amplitude distribution of WGMs; Bottom: GVD as a 
function of the resonant wavelength. For WGMs with transverse indices of (a) q = 1,2,3, 
p = 0 and (b) p = 0,1,2, q = 1. Dashed line: material GVDm curve. 


However, one can note that for the second radial order WGMs with q = 2, the GVD has been 
further reduced to 0.086 ps/km/nm, that is only 1/80 of the material GVDm and 1/35 of that 
for q = 1. In this case, Kerr comb generations can be significantly different, dependent on which 
radial order mode is excited. Comparing Fig. [4] with Fig. [3j one can see that the management 
of GVD by choosing WGMs with different q indices is more remarkable. Besides, we also find 
that the change of GVD due to the q index variation is much larger in the resonator with a 
smaller R. 

4. Experimental observation of significant different primary combs in a single MgF 2 
resonator 

To experimental verify that the geometry dispersion plays an important role in WGM res¬ 
onators, we have carried out experiments on a disk resonator made of MgF 2 . The resonator 
was fabricated by carefully polishing the rim of a disk preform with a radius of 6 mm and 
a thickness of 1 mm. It features an intrinsic Q factor of 10 9 . For a well polished resonator, 
WGMs with lower order mode indices can possess the same Q factors and cannot be distin¬ 
guished through Q factor measurements. It has been experimentally shown that one can choose 
the fiber taper to excite WGMs with the radial indices from q = \loq=\\ simultaneously (36). 
Hence, a fiber taper was used to couple the light into and out from the resonator as illustrated 
in Fig. [T] (b). The throughput of the fiber taper was then monitored with a photodetector and 
a high-resolution spectrum analyzer (APEX 2440B), respectively. The tunable cw pump laser 
with a sub-KHz spectral linewidth at 1552 nm was used. It was then thermally locked to ex¬ 
cite different WGMs (HD- Depending on the pump condition, both primary Kerr combs and 
full developed Kerr combs were observed. The value of single FSR spacing (vfsr) can thus be 
measured. 

The giant shift in dispersion predicted by our calculations should therefore allow for the 
possibility to generate primary combs with significantly different multiplicities, by pumping 
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Fig. 5. Experimental observation of three primary combs showing very different spacing 
of (a) 4-vfsr, (b) 164-vfsr and (c) 229 -vfsr- Three different WGMs were pumped in the 
same resonator, and the large difference of GVD explains the large difference of primary 
comb multiplicity. 


different radial modes in the same resonator. In particular, one should be able to observe very 
different primary comb spacing in the same resonator with a fixed incident pump power, es¬ 
pecially in the case of fiber taper coupling where many different radial order WGMs can be 
excited at the same time. Figure [5] present three experimentally obtained primary combs in the 
same MgF 2 resonator when three different WGMs were pumped. The spacing of 23.7 GHz , 
971.9 GHz and 1351.1 GHz corresponds to 4-vfsr, 164-vfsr and 229-vfsr, respectively. The 
corresponding ratios between the largest spacing and the smallest one is about 50. Although 
the detuning between the pump and the resonance can also change the spacing of the primary 
comb, the difference is usually limited due to the requirement of on-resonance pumping. In our 
case, we usually see a ratio of less than 2 when changing the detuning within one mode with 
a fixed incident pump power. Therefore, we believe that the geometry dispersion differences 
cause this observation, which is expected as previously discussed. 

5. Towards the mid-IR regime 

As shown in Table [I] BaF 2 has a larger zero material dispersion wavelength (Ao) of 1.93 gm, 
compared with other materials studied here. In addition, the GVDm of BaF 2 in the anomalous 
regime also grows much slower. This is an important signature indicating that it could be very 
favorable for the broadband comb generation in the mid-IR regime. As shown in Fig.[2j 2.5/im 
marks the first Kerr frequency comb generated in the mid-IR regime where a MgF 2 WGM 
resonator was utilized. It also indicates that comb generation is feasible for a material GVDm of 
25 ps/km/nm. It should be noted that BaF 2 has a nonlinear refractive index coefficient (ji 2 , in 
10 -16 cm 2 /W) of 2.85, which is larger than 1.90 of CaF 2 and 0.92 of MgF 2 (25 ]. If we assume 
a BaF 2 resonator with the same Q factor and the same geometry dispersion, one can thus expect 
a Kerr comb generation capability at A = 5.5/im. Considering that the demonstrated Kerr comb 
at 2.5/im in MgF 2 has a very high phase stability 0. the capability of generating such combs 
further into the wavelength above 5/im in BaF 2 is thus highly desirable. 

Here, we present the corresponding GVD calculation on WGMs made of BaF 2 for Kerr comb 
generations towards mid-IR. Figure [6] (a) shows the GVD for the fundamental WGMs in res¬ 
onators with different R and r = 0.1 mm. As BaF 2 has a material GVDm at 1.93 gm, it is thereby 
feasible to push the zero dispersion WGM resonance wavelength Ao close to 2.5 gm, which can 
potentially facilitate broadband comb generations into this regime. The corresponding Ao as a 
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Fig. 6. (a) Calculated geometry dispersion influence on the total GVD for a disk BaF 2 
with different R and r = 0.1 mm. Dashed line: material GVDm curve., (b) Calculated zero 
total dispersion wavelengths Ao for R from 0.5 mm to 5 mm. Inset: distribution of the 
electric field amplitude for a fundamental WGM with R = 1 mm and r = 0.1 mm. (c) 
Calculated radial mode order influence on the total GVD for a disk BaF 2 with R = 1 mm 
and r = 0.1 mm. WGM indices: q = 1,2,3,4,5,6 and p = 0. (d) The corresponding electric 
field amplitude distribution of WGMs. 


function of the main radii R are shown in Fig. [6] (b). The inset gives the field distribution of a 
fundamental mode at about 2.2 gm in a resonator with R = 1 mm and r = 0.1 mm. 

Figure [6] (c) presents the GVD as a function of WGM resonance wavelengths for different 
radial order modes in a BaF 2 resonator with R = 1 mm and r = 0.1 mm. The corresponding field 
distribution is shown in Fig. [6] (d). As one can see, GVD values at 3.5 gm for different radial 
WGMs are anomalous and reduced with the increase of q indices. In the case of prism coupling 
setup, one can change the coupling angle to selectively excite WGMs with high q values. The 
corresponding wavelengths and indices of those modes closest to the zero total dispersion in 
Fig. [6] (c) are listed in Table [2] For q = 1 and q = 2, the change in ZDW of about 240 nm is 
observed. ZDW beyond 3 gm is feasible in higher q order modes. 


Table 2. Indices of WGMs and zero total dispersion wavelength for MgF 2 with R = 1 mm 


and r = 0.1 mm 


Wavelength 
(Ao (WGM)) 

m 

q 

P 

Wavelength 
(Ao (WGM)) 

m 

q 

P 

2.1863 

4176 

l 

0 

2.8877 

3100 

4 

0 

2.4290 

3733 

2 

0 

3.1100 

2860 

5 

0 

2.6608 

3386 

3 

0 

3.3272 

2656 

6 
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6. Conclusion 


In conclusion, we have investigated the material dispersion profiles of different optical materials 
that support WGMs with Q factors above 10 9 . We have studied the effect of geometry dispersion 
on the total GVD of WGM disk-resonators. This analysis has enabled us to demonstrate that the 
major radii of WGM resonators and the radial order indices are key parameters for dispersion 
tailoring. In particular, we have shown that the dispersion shift induced by selecting different 
radial families of modes in the same resonator allows for the generation of versatile combs. 
Barium fluoride could benefit from this feature and it offers the best dispersion profile in the 
mid-IR spectral range. We expect this work to facilitate the efficient generation of Kerr optical 
frequency combs in the telecom and mid-IR regimes. 
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